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SUMMARY 

This Interim report consists of three separate reports: "Laser 
Altimeter Ranging Error versus Received Pulse Width and Amplitude 
for Leading Edge, Half Width, and Centroid Detection," "Final Report 
GAMES Master Oscillator," and "Sample Rate and Number of Bits 
Required for the GLAS Waveform Digitizer." 
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1. introduction 

Laser altimeters measure the time of flight of the laser pulses to 
determine the range of the target. The simplest altimeter receiver 
consists of a photodetector followed by a leading edge detector. A 
time interval unit (TIU) measures the time from the transmitted 
laser pulse to the leading edge of the received pulse as it crosses a 
preset threshold. However, the ranging error of this simple detection 
scheme depends on the received, pulse amplitude, pulse shape, and the 
threshold. In practice, the pulse shape and the amplitude are 
determined by the target target characteristics which has to be 
assumed unknown prior to the measurement. The ranging error can be 
improved if one also measures the pulse width and use the average of 
the leading and trailing edges (half pulse width) as the pulse arrival 
time. The ranging error becomes independent of the received pulse 
amplitude and the pulse width as long as the pulse shape is 
symmetric. The pulse width also gives the slope of the target. The 
ultimate detection scheme is to digitize the received waveform and 
calculate the centroid as the pulse arrival time. The centroid 
detection always gives unbiased measurement even for asymmetric 
pulses. In this report, we analyze the laser altimeter ranging errors 
for these three detection schemes using the Mars Orbital Laser 
Altimeter (MOLA) as an example. 

2. Leading Edge Detection 

The standard deviation of the ranging error for leading edge detection 
is given by 
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where at is the rms jitter in threshold crossing time and c is the 
speed of light. To obtain at , let us denote the photodetector output as 
ipd(t), and define the signal and the noise as 

^(0=<ipdU)> 

n(t) = ipd(t)-<ipd(t)> 


Let sth and Tr be the threshold level and the average leading edge 
threshold crossing time, approximate s(t) by the first two terms of 
its Taylor expansion about Tr. The threshold crossing time, //?, can be 
written as the solution to the equation 

Sth=s(tp) + n{tp) 

^s(Tp) + s\Tp ){tp-Tp) + n{tp ) ^ ^ 


Since s,h = s(Tp), (tp -Tp) ~ -n(tp)/s'(Tp), the variance of the 
threshold crossing time becomes 


crf(Ta)=<('R-T„f 




(4) 


with an^(t) the variance of the total noise at time t. 


Let us define the signal as the number of primary photoelectrons 
accumulated over the equivalent integration time of the lowpass 
filter, Tn- The signal can be written as^ 


s{t) = nGz„ jhiT)p(t-T)dr (5) 


with n the average number of detected photons per pulse, G the 
average gain of the photodetector (usually an APD), h(t) the impulse 
response of the lowpass filter before the discriminator, rh(t)dt = l, 


^ D. L Snyder, Random Point Processes, John Wiley & Sons, New York, 1975, p. 171. 
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and p(t) the normalized pulse shape function, j p{t)dt = \. The one 
sided noise bandwidth of the lowpass filter is given by 


^ -oo ^ 0 

The equivalent integration time of the lowpass filter is given by 

^ 2BW„ 

The derivative of the signal is given by 

5'(r) = nGT,jr^//(r)p'(^- r)r/T (8) 


For Gaussian pulse shape, 


P(0 = 



P'(0 


1 -t 

V2^cr^ <jI 



(9) 


where ctl is the rms pulse width. Here, we assumed the received 
pulses arrive at t=0. Note the rms pulse width is related to the FWHM 
pulse width, Tfwhm, by 


. _ T^FWHM _ q p g ^ FWHM 

^ 2V'21n(2) ■ 2 


( 10 ) 


The impulse response of the lowpass filter is determined by the 
actual filter used. The bandwidth of the APD and the preamplifier is 
assumed much wider than that of the lowpass filter. Appendix A 
shows how the impulse response function was obtained for a 5 pole 
Bessel lowpass filter. 

The threshold level, sth, is set such that the receiver false alarm rate 
does not exceed a certain level over a given range gate interval. 
Appendix B shows a detailed derivation of the optimal threshold level 
using MOLA as an example. The average threshold crossing time at the 
leading edge, Tr, is the first of the two roots to the equation 
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s{T^)-s„, = nGr„ jhit)piTf( - tMt-% =0. 

— <50 

The total noise at the discriminator input consists of three parts, the 
shot noise from the constant background illumination, the thermal 
noise of the preamplifier, and the shot noise due to the received 
signal. The one sided noise current spectral density at the input of the 
preamplifier can be written as2 

N„(/) = 2GVA* + W„„„ + 2GVi(r), (A^Hz) (12) 

where Af? is the constant photoelectron rate due to background 
illumination, Namp is the preamplifier noise current density, and F is 
the APD excess noise factor. The total noise variance can be written 
as 

(jlit) = jN„U)T„ ( 13 ) 

The standard deviation of the ranging error can be obtained by first 
calculating the threshold level, sth, as in Appendix B, solving (11) for 
the average threshold crossing time, Tpi, substituting it and (8) and 
(13) into (4), and then substituting (4) into (1). 

3. Half Pulse Width Detection 

The pulse arrival time in this case is defined as (tR+tF)/2 with tR and 
tp the threshold crossing times at the leading and trailing edges. The 
standard deviation of the ranging error in this case is given by 

_2 ^(Tg)+af(TF) _i I (7^(Tg) a;(Tr) ,,,, 

where Tr and Tf are the average threshold crossing times at the 
leading and trailing edges and they can be obtained by solving (11). 
Compared with (1), the rms jitter for half pulse width detection is 

2 R. G. Smith and S. D. Personick, "Receiver design for optical fiber communication 
systems," in Semiconductor Devices for Optical Communication, Spring-Verlag, New 
York, 1980, ch. 4. 
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improved by 1 /V 2 over the lead edge detection when the jitters at the 
rising and trailing edges are equal. The noise variances, OniTpi) and 
<yn(T^R ) . and slopes, s'(Tr) and s'(Tf) can be obtained by substituting 
Tr and Tp into (13) and (8). 

4. Pulse Centroid Detection 


The detected pulse arrival time using centroid detection is the same 
as that of the half pulse width detection provided the received pulse 
shapes are symmetric. The standard deviation of the ranging error for 
centroid detection is given by^ 


- 2 , 17’ 


2 f 


FXJ- + 


V 


2GW 


(15) 


where T is the observation time. Nth {A^/Hz) is the one sided input 
noise current spectral density of the preamplifier, and q is the 
electron charge. We used T=Tfwhm in the subsequent numerical 
calculation. In theory, there should be an optimal value for T which 
minimized the effect of the background and thermal noise and best 
preserve the received pulse shape. We will study this issue in greater 
detail in the near future. 


5. Numerical Calculations and Discussions 

Appendix C shows a Mathcad'l’M program for numerical calculation of 
the ranging error versus received pulse width and amplitude using 
MOLA parameter values as an example. The results are presented in 
Figures 1 through 4. 

Figure 1 shows a plot of measurement bias as a function of the 
received pulse energy (amplitude square) for a 60ns FWHM pulse width 
and a 5.54 MHz lowpass filter bandwidth. The measurement bias of the 
leading edge detection decreases as the pulse energy increases. In 
theory, this bias may be corrected if one also monitors the received 
pulse energy. The measurement bias for half pulse width detection is 

3 C. S. Gardner, "Ranging performance of satellite laser altimeters," IEEE Trans. 
Geoscience and Remote Sensing, Vol. 30, No. 5, pp. 1061-1072, Sept. 1992 
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almost constant, which equals to the lowpass filter propagation delay 
and can always be calibrated out. Note we have assumed the received 
pulse shape is Gaussian and centered at t=0. Because the impulse 
response of the 5 pole lowpass filter is not perfectly symmetric, 
there is still a slight increase in the bias over pulse energy for half 
pulse width detection as shown in Figure 1. 

Figure 2 shows a plot of measurement bias vs. received FWHM pulse 
width under constant received pulse energy and the same lowpass 
filter bandwidth. The half pulse width detection again gives a nearly 
constant measurement bias, which can be calibrated out once the 
lowpass filter propagation delay is known. On the other hand, the 
measurement bias vs. pulse width for leading edge detection first 
decreases because the pulse gets wider and then increases because 
the peak of the pulse approaches the threshold level. The rms jitter 
(error bars) also increases with the pulse width for both of the 
detection schemes. 

Figure 3 shows the standard deviation of the ranging error in 
nanoseconds vs. the received pulse energy for leading edge, half pulse 
width, and centroid detections. As expected, the half pulse width 
detection gives a lower rms jitter than the leading edge detection and 
the centroid detection always gives the smallest jitter. However, the 
rate of decrease of the rms jitter becomes slower and slower for the 
leading edge and half pulse width detections because the slopes of the 
signal at threshold crossings become shallower and shallower under a 
constant threshold level. 

Figure 4 shows the standard deviation of the ranging error as a 
function the received pulse width under constant pulse energy. The 
rms jitters increase very rapidly for relatively large pulse width as 
the threshold gets close to the peak pulse amplitude where the slope 
is shallow. There is a moderate improvement in the standard deviation 
of the ranging error by using centroid detection. 
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Appendix A. impulse Responses of Fifth Order Bessel Lowpass Filters 

Xiaoli Sun/JHU, May 9, 1 995 


Pole locations: (De Verl S. Humphergs, The Analysis, Design, and Synthesis of 
Electrical Fitters, Prentice-Hall, 1970, p. 421) 


k ;= 1 .. 5 

s, :=-0.9576766+ 1.4711244j Sj :=-0.9576766 - 1.4711244j 
S 3 := - 1.3808774 + 0.7179096j :=- 1.3808774 - 0.7179096J 

S5 :=- 1.5023160 


Multiplication factor (such that H0)=0) 

-1 

aO = 11.2128369 


Coefficients of partial fractions 



si -= Sj 


' ( si — Sj )■( Sl — S3 )•( si — S4 )•( si — S 5 ) 

si '= $2 


Aj = 0.6193646 + 1.0986249J 


Aj •= 


^ 

( si — Sj )•( si — S 3 )■( si — S^)-( si — Sj ) 


A 2 = 0.6193646- 1.0986249J 


si Si 


A3 := 


aO 


( si — Sj )■( si — Sj )•( si — S 4 > ( si — Sj ) 
si :=s^ 

aO 


A4-= 


(sl — Sj)*(sl — S 2 )‘(sl — S 3 )*(sl — S 5 ) 


si :=S5 


A, := 

5 (Sl-Sj)-(SI 


aO 

Sj )•( si — S3 ) (S 1 — S4) 


A 3 = -4.916811 -2.6131065j 


A 4 = -4.916811 +2.6131065j 


A 5 = 8.5948928 



The Impulse Response Function 


hh( t , BW “ 


exp( Sj^ ( 2 It BW 3^3 ) t) 
k 


2 n BW 


3dB 


The real impulse response function: 

h(t,BW3jg) :=Re(hh(t,BW3jjg)) 


Note the imaginary part 
is not canceled out completely 
due to computer quantization error. 


Numerical results of the normalized filter 
BW 3J3 :=^ -L (1 rad/sec) 

k -=0.. 500 tO^ -=0.02-k-sec 

hO^ :=h(t0j^,BW3^) 



0 2 4 6 8 10 

FWHM pulsewidth; ‘°k 

tl :=l sec Initial guess 

T1 ;= root( linterp( tO , hO , tl ) — 0.5- max( hO ) , tl ) 
t2 •= 3 sec 

T2 r=root( linterpf t0,h0,t2) — 0.5- max( hO ) . t2 ) 


WO := T2 - Tl WO = 2.204 *sec 

The actual pulse width: 

e.g., for MOLA 20 ns filter (3dB bandwidth=16.6MHz). 

T p\yjj\{( 16. 6‘ 10^ ) = 2.1 13* 10 "sec 

One sided noise bandwidth: 

j I* 100 - sec 

BW„:=-- h(t,BW3^)2(it 

• 0* sec BW 


= 1.039 




Appendix B. Calculation of MOLA 

Discriminator Threshold Level 

Xiaoli Sun/JHU, May 9, 1995 
Average # of background radiation photons/s 
I solar 1 Solar irradiance at the target (W/m^2/nm) 

1=0.26 Target diffusion coefficient 

^ tel 0-508 Receiver Telescope Diameter (m) 


^ obs •” -19 


Telescope obscuration ratio 


0 FOV 0.847- 10 
^ rcvr ■“ 0-63 

^ atmo •“ 0.5 
APD '“0.3 
AA. :=2.2 
R :=410- 10^ 

1 b := 50 lo" 


Receiver FOV (rad) 

Receiver optics transmission 
Atmosphere Transmission Coefficient (one way) 
APD quantum efficiency 
Receiver optical bandwidth (nm) 

Orbit altitude (m) 

APD bulk leakage current (A) 


Background noise power 


^ b * solar 


. ^ tel . 

"’^”1 2 ~) ~ 


■ ( I 'n obs ) 
2 


A. h •=■ 


^ APD Pb 


1.171.610"^^ 1.6- lO"^^ 


APD Spectral Noise Density 

G := 120 Average APD gain 

k := 0.007 APD ionization coefficient ratio 


P'=G-lCefiF+ (1-k^) APD excess noise factor 

NaPD-=2 g2 f A.b 


Pj, = 3.3* 10“^ 

A. b = 5.6' 10^ 

8 

= A.IMO** 

c.,n-19 


F = 2.82 


APD noise current density (A'^2/Hz) normalized wrt q 



Preamplifier Thermal Noise Spectral Density: 




195- 10“^) 


Ny : = 

24- lO'^ 
R , := 20000 


Preamp dark noise density given by EG&G (V'^2/Hz) 
APD load resistance (ohm) 


Equivalent noise current at input (A/Hz^1/2) 




= 1.99*10 


N. 


N 


amp 


R (l.6 lO'^^) 


2 Equivalent noise current density normalized wrt q 


Noise due to APD surface leakage current: 


I/.= 1510 


,-9 


APD surface leakage current (A) Nj^ :=2 


1.6- 10 


-19 


Other system parameter values: 

Range gate (km) 


Rg .-20 


T :=2 


R Q lOOO 


3- 10' 


8 


Range gate interval (sec) 


T = 1.33*10 


Matched Filter Noise Bandwidth (5 pole Bessel LPF) 


BW 3^ :=5.54- 10^ 

(Channel 2) 

BWj, := 1.039- BW3 jJ3 

(see Appendix A) 

1 

2-BW„ 

equivalent integration time for noise photons 

Equivalent average # of noise photoelectrons: 



F0= , 

0 = 4-87' 10"^ 

G^-F 



X = 8.69* 10 * 


The variance of the preamplifier noise and APD surface leakage current noise (Gaussian): 


a 



a = 2.59*10^ 


— = 2.16*10^ 
G 



The p.d.f. of the APD output: 


4 


Soo:=^g2f-ho Sqo = 4.44*10^ 


P pDo( z) ■= 


1 

1 



-z2 


IT, 



2- 

G(F- 1) z' 



r o(F-ii^i 

|2 

. 

*00 

. 


00 


(see the attached notes for the derivation of the above equation) 


Probabilities of false alarm: 


s 


0 •“ 



s 0 = 5.14- 10^ 


Standard deviation of the total noise 


z 


0 •“ 


"G p 0 
* 00 


Zq = -1.31’ 10^ 
z j :=30 


lower limit of the integral 
Upper limit of the integral 


Pfe(n-p) 


P pqq( z) cnorm 


|s 00 z-nj SQj 


dz 


(see the attached notes 
for derivation) 


The built-in function ■cnorm(u)' is the culmulative probability 
distribution function of a Gaussian r.v. with zero mean and unity variance 


Pfa^"T^ 1 - exp(-T BW 3 dB P ^(n j)) 


Assuming the average pulsewidth is 1/BW. 



M ;= 10 


k :=0, 1.. M 

NTk :=3 + 0.2 k Normalized threshold level 
TOL := 10 The integration error tolerance 

pfek •= P NTk ^ 

PFAk •= 1 — exp( "T- BW pfa^ ) 



Normalized Threshold level at PFA=0.01: 

PPFAk:=PFAM_k NNTk:=NTM_k 


Thrc 1= linterp( PPFA, NNT, 0.01 ) 

Thre = 4.5 
Sq = 5.14* 10^ 



for 




Notes to Appendix B: 

The average rate of false alarms can be calculated as^ 

'T'— 1 _ Pfa 

‘FA-x^ 


(1) 


where pfa is the false alarm probability at a fixed time and tw is the 
average noise pulse width above the threshold. It may be assumed that 
xw=1/BW3cib- The number of false alarms over a given range gate time, 
Tg, follows a Poisson distribution and the false alarm probability is 
given by 


PpA=l- exp|- ^]=1- exp[- ToBWjdePfJ 

The false alarm probability at a fixed time can be written as 


Pfa=| p(yi^o)dy= 

JYt 




(2) 


(3) 


where p(ylpo) is the p.d.f. of the APD preamplifier output and is equal 
to the convolution of the p.d.f. of the APD output and the p.d.f. of the 
preamplifier noise. The latter is assumed to be additive and Gaussian 
with mean, pg, and standard deviation, a. Usually, the preamplifier 
itself has zero DC offset, pg=lst/q with Is the APD surface leakage 
current, x the equivalent integration time of the lowpass filter, and q 
the electron charge. The variance is given by 

aNN3^+2I^BW„x2/q2 (4) 

with Namp the one sided input noise spectral density of the 
preamplifier in (A^/Hz) and BWn the noise bandwidth of the lowpass 
filter. 


The p.d.f. of the APD output in number of photoelectrons is given by2 

- (x-Gpo)^ 


PpoCxIPo)- 


1 


2jtG^F^o 


1+ 


(x^^o)(F-l) 


GF^o 






(5) 


^ M. I. Skolnik, Introduction to Radar Systems, McGraw-Hill, New York, 1962, p. 31. 
2p. p. Webb, R. J. McIntyre, and J. Conradi, 'Property of avalanche photodiodes," RCA 
Review, Vol. 35, pp. 234-278, June 1974. 



where \iq is the average number of primary photoelectrons given by 


,, _T1aPdP 
fio-— hf— +nq- 


( 6 ) 


with tiapd the APD quantum efficiency, Pb the received background 
noise power, hf the photon energy, and lb the APD bulk leakage current. 

interchange the order of integration in Eq. (3) and substitute 
u=(x-y)/o, and z=(x-Gpo)/(G2F^o)^^^ . 


Pfa= 


</2kc 


(y-x-fi 


^yj 


£1. 


- ~~ 2q2 fy >PpD(xlPo)dx 


I L Tlx ^]du|ppD(xlpo)dx 


Vo^Fllo 


PpDz(zlFo)-y|^ 


1 


1 + 


G(F-l)z 


\/g¥ho 

A normalized threshold is defined as 


i3^*P 




-z-‘ 


G(F-l)z 

>/G"Fpo 


Y y-r-(Gpo-bPg) 


(7) 

(8) 
(9) 


( 10 ) 


that is, the actual threshold is Yj times the total noise standard 
deviation plus the average DC component of the APD output. Usually 
Pg«Gpo, we can drop pg in the calculations. 



Appendix C. Calculation of Ranging Errors of MOLA 
vs. Pulse Width and Amplitude 
for Leading Edge, Half Width, and Centroid Detections 

Xiaoli Sun/Code 924/JHU, May 11, 1995 


System parameters 

Received laser pulse energy 

APD quantum efficiency 

Average # photoelectrons/pulse 

Received pulse width (FWHM) 

Matched Filter 3dB BW 
Average APD gain 
APD ionization coefficient ratio 
APD surface leakage current 
APD bulk leakage current 

Preamp input noise current density 


Ej.-=.51010 joule 
Ti := .30 

"bar'“ -JO 

(1.17- 1.6- 10 ^)- joule 

^FWHM '=60 10 ^ sec 

BW 5. 54' 10^' Hz 

G ;= 120 
k^:= 0.007 
[5 ;= 15 10~^ amp 
Ij, :=50 lO’^^ amp 



“bar = 817 



sec 


Backgound illumination power into the APD P b ^o”^‘ watt 
Noise and threshold levels (copied from Appendix B for Ch. 1): 

s 0 •- 5.14- 10 Standard deviation of the APD and preamp noises 

Thre .= 4.50 Normalized threshold for PFA=0.01 per range gate daytime 


1 

2-1.039 BWppp 


Equivalent noise integration time of the lowpass filter 


Background photoelectron rate: 


A. h •=• 


^ lb 

1.17- 1.6- 10 joule 1.6- 10 ^^ coul 


Total additive Gaussian noise density; 
Laser pulse shape function: 


2 1 


Namp-=Na + 


1.610"^®-coul 


RMS pulse width 
Laser pulse shape 


® l(T FWHM ^ 


^FWHM 

2-^2- ln(2) o L(T FWHM ^ = 2-548* 10~* ‘sec 


® L< T FWHM > I ^ ®L(TfWHM)^, 



Impulse Response of the lowpass filters: 

5 th order Bessel LPF 
Pole locations k:=i ..5 

Sj :=-0.9576766 + 1.4711244j Sj :=“0.9576766 — l.4711244j 
S 3 :=-1.3808774+ 0.7179096J -=-1.3808774 - 0.7 179096J 

S 5 :=- 1.5023 160 


Coefficients of partial fractions 

Aj := 0.6193646+ 1.0986249j A 2 := 0.6193646— 1.0986249j 

A 3 ■=-4.916811 - 2.6131065J A^ := -4.916811 + 2.6131065j 

Aj := 8.5948928 

Impulse Response: 


^A|j^ exp( s,^ ( 2 Jt BW Lpp) t) 

* 2* 7C* fifW j^pp 

k 



Average signal (secondary photoelectrons p>er x second): 


s( t, n bar T fwhm ^ 


' “ bar ^ n 


-5 




^ FWHM ) “ u ) du 


Total noise variance calculation: 

APD excess noise factor F :=kg^G+ ^2 — — j •( 1 — k 

variO.nbar.TpwHM) n bar"' n^‘ f P{ T. T ) h( t - x dx 

° bar' FWHM ) 1^ bar FWHM ^ ® 0^ 


Derivative of the signal: 


*^l(Tpwhm) 


dsdKt.nbar.Tp^^^jjjj^) :=nbar Gx 


r P^*'^FWHM^ 


-5<TL(Tp^fU^() 


“ » T FWHM t " ** ) du 



Example (ATU signal laser): 

1 -=0..38 tj •=(i) 0.05-(Xn+Tpyvjjj^)- l O T 


Si :=s(ti,nbar.T 


FWHM 


) 


THRE- *= Thre* s q 



Solve for the threshold crossing times 

equation to solve H T n tl ) := root( s( tl- sec, n T )- Thre s q, tl ) 


Rising edge: 


1 


tl — initial guess 

TR =KTp;^^jjj^,nj^,tl) 

Falling edge: 

t2 .= (0.5 Tp^yjjj^+ ^ n) — 

^ ^ FWHM’ “ b»’ ) 

Half width location: 

. TF+TR 

TO .= 

2 

RMS jitter at threshold crossing: 


TR = 1.03*10 


.“8 


TF = 1.303*10 


TO = 7.029* 10 


® xr •“ 


var( TR sec, □ T pyyjjj^ ) 


dsdt( TR sec. n bar ^ FWHM ) 
xh -0.5- J~t 


I var(TF sec,D(^,Tp^yjj^) 

Oxf- 1 

dsdt( TF- sec, n T ) 


xr^ + o 


•= (- l( T fvvHM ~ 

“bar n 


tar^ ’’ I 2 ay 


^ xr ” 4*836* 10 *sec 


o = 3.507* 10 *sec 


^ xc ^ 2.61* 10 *sec 




NUMERICAL CALCULATION -- Ranging error vs. pulse energy: 

Received Photons/pulse: 

' ■“0 .10 nbarj :=( i- 0.25 + 0.25)- n (jgj. 

Solve for the threshold crossing times 
equation to solve 


FWHM’ " bar tl sec. n T p^HM ^ ~ s 0’ ^ 


Rising edge: tl. :=- 


FWHM+ ^ 

L 


-T 


■('-3 


•0.2- — initial guess 
sec 


Falling edge: t2 := ( 0.5- T pwHM ^ n) ~ inifisl guess 

sec 


TFj •=f(Tpyyjj^,nbarj,t2) 

Average (half width): 


THj :=0.5 (TR. + TF.) 

RMS jitter at threshold crossing: 

var( TRj sec.nbarj.Tpvv^) 

OTj .= • 

y dsdt( TRj- sec , nbar. , T ) 

:=-i-|(orp2 + (<jf.)2 


var( TFj- sec, nbatj , T pwHM ^ 

, **(TFj Sec,nbaTj,Tp^rHM^^ 


acj .= 


nbaf; 


•Ol(T 
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NUMERICAL CALCULATION -- Ranging error vs. pulse width: 


Received pulse width: i -q 14 

Solve for the threshold crossing times 
equation to solve 


TFWHMj •— (i'20'10 ^ + 60' 10 ^)'sec 


^ FWHM’ " bar ) •“ tl- sec. n T pwHM ^ * 0’ ^ 


Rising edge: 


tl- :=(-0.5 TFWHM- + T „)■ — 

sec 


initial guess 


TRj :=f(TFWHMi,nbar,tlj) 

Falling edge: t2j :=(0.5 TFWHM. + t „)■— initial guess 

sec 

TF. :=f(TFWHMi,nt^,t2i) 

Average (half width): TH- := 0.5- ( TRj + TFj ) 

RMS jitter at threshold crossing: 

var(TRj sec.nj^.TFWHM-) I var( TF- sec, n - , TFWHM ) 

OTj = — of. := — ^ 

J dsdt(TR. sec.njjg^.TFWHM;)^ J dsdt( TF- sec, n TFWHM^ 

<riij J (aTj)^ + 
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1 Introduction 

The Gravity and Magnetic Earth Survey (GAMES) is a mission to measure the 
Earth’s gravity field by recording the relative velocity between a main satellite and 
a sub-satellite which are in identical orbits. The relative velocity wiU be measured 
to about 20 fim/s when the range is 1(X) km and 50 fim/s when the range is 200 km. 

There will be three optical ranging signals. A pseudo-random binary sequence 
at 3 1 .25 Mbits/s . will intensity modulate the tracking laser diodes. A mirror retro- 
reflector on the sub-satellite returns laser light from the main satellite back to 
it Correlation of the received light with the transmitted signal will determine the 
inter-satellite distance to 4.8 meters. A second laser will be intensity modulated by 
medium and fine ranging sine waves. The medium signal will also be 31.25 MHz. 
Resolution of the difference between the transmitted and received phase to 5 de- 
grees will determine the inter-satellite distance to 70 mm with a 4.8 meter range 
ambiguity. The fine ranging signal, 2 GHz, has a 75 mm range ambiguity. As- 
suming the ranging frequency is known exactly, determining the transmitted and 
received phase difference to 0.1 degrees will yield the range to 20 /zm. The phase 
difference will be reported about 250 times each second. Tracking the phase dif- 
ference as a function of time will give the inter-satellite relative velocity. Velocity 
variations are expected to take longer than 10 seconds. 

This report describes the system oscillator specification required to meet the 
performance goals. 
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( 1 ) 


The transmitted fine and medium ranging signals may be written as 

i;7(/) = uo sin(f(Wc + 0o) 

where o)e is the modulation angular frequency and ^ is the initial phase. The 
received signal is an attenuated, time shifted copy of the transmitted signal and can 
be written as 

v/f (t ) = Vq - Trf)n>c + 4>o) (2) 

where X 4 is the round trip time delay. Denote the transmitted and received phases 
by 


arfO = tcoc + <l)oit) (3a) 

ccrU) = it-Td)(Oc + 4>oit-Ta) (3b) 

where account for oscillator imperfections have been taken by writing a time 
dependent phase. Multiplying the transmitted and received signals yields a term 
proportional to the sine of the phase difference 

Aa(r) = arit) - a/t(0 = XdO)c + ^(0 - M - tj) (4) 

The measured phase difference is proportional to the time delay and firequency and 
also depends on the oscillator phase variations. 

Because the phase difference is proportional to the riequency, the accuracy 
of the fine range measurement depends critically on the stability of the master 
oscillator. This report summaries the oscillator frequency and phase variations 
effects on the fine ranging performance. The analysis focuses on the fine ranging 
system, but the same mathematics governs the medium ranging system. 

Measurement of the phase difference, Aa(r), could be p^onned by mixing 
the transmitted signal with the received signal. While doing so is functionally 
simple it has the disadvantages that the final output signal is essentially DC and 
subject to 1// noise, DC drifts, and other measurement errors. In fact the received 
signal is mixed with a local oscillator offset from the transmitted signal frequency 
by an intermediate frequency, o),-. The resultant low frequency signal is digitized 
so that the final phase determination is immune to 1// noise, DC drifts, offset 
errors, etc. 

The rest of this work is structured as follows. Section 2 summarizes the 
GAMES oscillator performance requirements. Section 3 describes crystal oscilla- 
tors and oscillator frequency and phase noise. Methods for expressing oscillator 
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noise arc covered. Also covered is the phase noise power spectrum after taking 
the difference of oscillator phase taken at different times. Section 4 relates the 
specification of a few commercially available oscillators to the requirements. The 
appendix derives the expressions used to determine the oscillator performance 
requirements. Finally there is a short list of useful references. 


2 GAMES Oscillator Requirements. 

There are three constraints on the oscillator performance. The desired relative 
velocity resolution limits the permissible rate of frequency change, both determin- 
istic and random. Oscillator phase noise introduces random noise somewhat like 
amplifier thermal noise. The permitted noise level must be specified. The absolute 
frequency stability determines the inter-satellite range accuracy. 

For GAMES, the relative velocity between the main and sub-satellite must be 
measured to 20-50/im/s. So that the oscillator will be an insignificant error source, 
it will be specified to introduce a velocity error no larger than 1/i.m/s and phase 
noise no greater than the thermal noise contribution of the PMT pre-amplifier. 
Table 1 defines the symbols used in this work. 

The required oscillator specifications may be summarized by 

1. A range rate accuracy of Au measured over a time period Ar requires two 
range measurements with a permitted range error of Ar = At Av/2. The 
maximum permitted frequency change is given by A/c//e < A/Au/2r 
or A/c/(/cAr) < Av/2r = 2 x 10~‘Vs, 1 x 10”*/hour, 2 x lO'Vday, 
7 X 10“Vmonth, or 7 x 10"^/year. 

2. So that the phase noise power spectral density to carrier ratio is about the 
same as the amplifier thermal noise power spectral density to signal ratio, 
the 2 GHz phase noise for offset hnequencies greater than 1 kHz must be less 
than -65 dBc/Hz. This specification also fixes the low frequency noise as 
will be shown. 

3. A frequency accuracy of Mdfc yields a range accuracy of Ar/r. Achiev- 
able frequency stabilities of 10~^ over the mission life of GAMES limits the 
absolute range accuracy to about 2 mm over the mission life time. 

The range rate requirement is challenging, but should be possible. Typical long 
term drift rates are about lO'Vyear. Sami Asmar, JPL, says frequency stability 
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fc = 2 GHz 

exactly 

(Oc = 2jtfc 

exactly 

fc + A/c 

be the actual oscillator frequency 

Au>c = 2jtA/c 


r(0 

be the actual range 

v{t) 

be the actual velocity difference between the crafts 

Ar 

be the required absolute range resolution 

Av 

be the required absolute velocity resolution 

ffr(t) 

transmitted phase 

a«(0 

received phase 

Aa 

difference between transmitted and received phase 

Aa, 

change in Aa due to a range change 

Aa/ 

change in Aa due to a frequency change 

c 

The speed of light in vacuum. 


Table 1: Symbols used in this section. 


of 10“*^ may be possible over several tens of seconds. Reference [2], page 219 
repOTts aging rates of 10~* Vday. 

The phase noise specification is satisfied by most oscillators, even after fre- 
quency multiplication to 2 GHz. Note that this requirement implicitly specifies the 
low frequency noise as is discussed later. 

The absolute frequency stability is a goal, radier than a requirement It may be 
possible to measure the oscillator frequency with the GPS system and correct for 
long term drift. Absolute ranging would then be improved if required. 

See appendix for derivation of these requirements. 

GAMES will have a calibration path which in principle could be used to detect 
frequency changes. Substituting the calibration path time delay, Tc, for die range 
path delay 2r/c in equation (24) gives the calibration phase shift as 


Oc + Atte = 2Tl[fc -I- A/c]Tc 


( 5 ) 


The change in calibration path phase due to a frequency change is then given by 
A«c = 27T A/cTc. The ratio of calibration phase change to range phase change is 


Aac _ 
Att/ ~ Tj 


( 6 ) 
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where Xd is the inter satellite round trip range time delay. The calibration time delay 
is likely to be 10 ns to 100 ns. The range delay is about 1 ms. Therefore the cali- 
bration channel will then be 10~^ to 10"* times less sensitive to frequency changes 
than the ranging channel and using the reference path to correct for frequency 
variations is not feasible. 


3 Oscillators 

3.1 Introduction to Quartz Crystal Oscillators 

Quartz crystal oscillators are used as the frequency selective filter in oscillator 
circuits. Because quartz is piezo-electric, an oscillating electric field can generate 
acoustic waves which can in turn generate an oscillating electric field. The res- 
onator Q can be very large and is highly frequency selective. Quartz crystals can 
operate over the frequency range 100 kHz to 150 MHz.[l 1, Page 57] The crystal 
may be cut in one of several ways, each with advantages and disadvantages. The 
most common cut is AT, in which the crystal is cut once along the crystal axis. 
AT cut crystals are inexpensive but not as stable as SC (stress compensated) cut 
crystals which require two cuts in directions not parallel to any crystal axis. SC 
cut crystals have reduced temperature and vibration sensitivity, better short-term 
stability, lower aging rates, and higher immunity to radiation than AT crystals[2. 
Page 227]. 

As a quartz crystal is used, its frequency changes with time. There are three 
major components to the change, called drift or aging[2, (Chapter 12]. After being 
off, there is a frequency drift as the crystal reaches thermal equilibrium. This 
takes a few tens of minutes to hours. Then over the next few days, the frequency 
changes as impurities, which condensed on the crystal while it was off, are driven 
off. Hnally, the frequency drifts as the crystal relaxes. When the crystal is new, 
surface stress induced by crystal cutting relax causing ftequency changes. Also 
the internal inter-atomic restoring forces weakened with use. This drift can be well 
characterized. It tends to be rapid when the crystal is new. After one or two years, 
the aging rate decreases to a steady, low value. In GAMES, only the last drift 
effect will be important since the oscillator will be run continuously. Drift rates 
and high ftequency phase noise can be traded for each other. See the paragraph 
before equation (12) 

Generally low frequency crystal oscillators have long teim drift rates (aging 
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rates) lower than high frequency oscillators because the lower fiequency crystals 
tend to be larger. With more material, the restoring forces are distributed over 
more atoms and the weakened bonds have a proportionally smaller effect Large 
crystals can be used as high frequency resonator by operating on a vibration mode 
above the fundamental, known as overtones. Therefore many crystals are operated 
on the third or fifth overtone. It should be noted that overtone frequencies are not 
harmonics of the fundamental 

3.2 Oscillator Noise 

In addition to long term frequency drift, oscillators have random fluctuations in 
frequency, phase, and amplitude. Here amplitude noise will be ignored. Let the 
output of an oscillator with angular frequency be described by either 

v(t) = uosin(J cDoit'W + <l>o) (7a) 

v(0 = vosin(o)ot + 4>(0) (7b) 

where uo is the amplitude. Equation (7a) describes the oscillator imperfections in 
terms of fluctuations in frequency, ru«(r). Equation (7b) describes die oscillator 
output is terms of an average frequency, a}„, and random phase fluctuations, 0(r). 
Generally (7a) is used for low frequency variations, less than 0.1 Hz, while (7b) 
is used for high frequency variations, greater than 100 Hz. Either may be used in 
the transition region. 

Oscillator frequency noise may be described in terms of the Allan variance. 
Let the oscillator fi^uency be measured during each of N adjacent time intervals 
of length T. Let /* be the frequency measured during the k’th interval. The Allan 
variance is (tefined by 


Hie Allan variance can be normalized u> the carrier frequency by dividing by the 
square of the average frequency. The Allan variance can be converted to frequency 
space and plotted in terms of frequency offset from the carrier. Generally, however, 
it is left as a function of the measurement intervaL r, which typically ranges from 
a few tens of milliseconds to a thousand seconds. Aging frequency changes are 
removed before the Allan variance is computed. 
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Oscillator noise may instead be described in terms of phase noise by starting 
with (7b). The sin( ) term is expanded to yield 

u(r) = uo {sin(£D^,r) cos{(f> (t)) + cosi<Oot) sin(0 (r))} (9) 

Assuming that <p{t) is always small, cos(^ (r)) can be approximated as unity while 
sin(^(r)) can be approximated by its argument Thus (9) can be written as 

vit) « i;o sin(a>o/) + uocos(a)<,r)0(r) (10) 


The power spectrum of u(t) can be computed as 

. vlhd{(Oo + 0i) . vlSd((Oo-(o) . \^S^io>o±(o) 
V{(o) « ; 4- T + r 


where 8d (•) is the Dirac delta function and (to) is the phase noise power spectrum. 
Note that this is an approximation that depends on 0(r) 1 and therefore can not 

be used at frequencies very near the carrier. 

Typically the phase noise is specified as the ratio (in dB) of noise power in a 
bandwidth to the carrier power, Vq/4 at a frequency from the carrier, a>„=(o±Q>o. 
This is known as dBc. 

Conversions can be made between phase noise power spectral (knsity and 
frequency noise power spectral density as described in [1 1, C3iapter 9]. 

The phase noise power spectrum can be divided into three regions[l 1, Section 
5.3) as shown in Figure 1 . At ftequencies far from the carrier, large o>„, the power 
spectrum is flat due to flat, thermal noise in the oscillato- amplifiers. Generally this 
noise to carrier ratio can be decreased by running the oscillator at higher power 
levels. However, higher power levels in the crystal tend to increase the crystal 
aging rate because of the larger vibration amplitudes. At fircquencies closer to the 
carrier, the phase noise to carrier ratio is given by 



where o)m is the fiequency offset from the carrier and Q is the crystal resonator 
Q[1 1, Page 53). At frequencies even closer to the carrier, the phase noise increases 
as l/o)^duetoflickernoise,[ll,Page53]randomwalkoffrequency. The Transition 
frequency, the intercept between the 1/w^ and lines, depends on the flicker 
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Figure 1: Typical phase noise power spectnun showing ihe three regions: 1//^, 1//^, 
and the flat, white noise pans. 


noise of the circuit components. It should be noted that sometimes the white or 
flicker parts of the spectrum are so large that the 1//* part does not appear. 

The power spectral density in each region can be written as 


S^i(Om) = fl. 



(13) 


where n = 0, 2, or 3 and is the power spectral density at frequency a>m = <Um,- 
Suppose that <o,^ is the angular frequency at which the flat spectrum intercepts the 
1//^ spectrum (or the 1//^ if there is no 1//^). In Figure 1, q>^ is 2n 10 kHz. 
Then the spectrum can be bounded by 


S^{(o„) < 


Oo ,<Om>(Omo 

«0 (2^) . < <0«o 


(14) 


This proves the surprising result that specifying the high fireqiwncy phase noise 
bounds the low frequency phase noise. 

Recall that the entire phase noise analysis depends on the assumption of small 
amplitude phase noise. At some low frequency, often it is in the 1 Hz to 10 Hz range, 
this assumption breaks down and the system is no longer linear as is required for 
a Fourier description of the system. Low frequency phase noise is best described 
in the time domain. 
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Frequency multiplication by N increases the phase noise power by N^, the 
same as increasing the oscillator frequency. Therefore, assuming the crystal Q is 
frequency independent running the oscillator at a higher frequency has the same 
impact on phase noise as multiplying the frequency. In general, however, the Q is 
frequency dependent, decreasing at high frequency, and there may be an advantage 
in frequency multiplying. 


3.3 Power Spectrum of Phase Differences 


As described in the introduction, the distance measurement is computed from the 
phase difference Aa(r) which depends on ^(r) and — Xd). Qearly the two 
phases are correlated so the power spectrum of the phase difference can not be the 
same as the phase noise itself. It can be shown[ll. Section 10.1] that the phase 
difference power spectrum, 5^ is related to the phase noise power spectrum, 
S^Uom), by 

= 4S^(o>„) sin^ (15) 

It is convenient to bound the power spectral density by using the relations that for 
aU X, sin^x < x^ and sin^x < 1 which yields a bound given by 


Si^aicOm) < 


• l^ml — 

4S^(to) , > 2 tj ‘ 


(16) 


Substituting (13) into (16) yields 


SsaiO)^) < 


4S^(tu) 


. I<U„| < 2xj ‘ 
, |<u„| > 2 t7‘ 


(17) 


where, as before, /i = 0, 2, or 3, is the exponent appropriate for the frequency. For 
frequencies less than the inverse delay time, the phase difference power spectrum 
is r^uced by differencing. Unfortunately, the low frequency, n = 3, spectrum 
still grows as l/co„. As will be shown, this effect is not important for velocity 
measurements although it is important for range measurement 

The phase difference variance, is bounded by the integral of (17) over all 
frequencies. 

cfL<2rS^aMda>^ (18) 

Jo 

Unfortunately, the integral in (18) is infinite. As usual, a high frequency cutoff is 
imposed by the system electronics. But the 1 /o)^ factor still generates a divergence 
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at low frequencies without an obvious cutoff frequency except ftM* the breakdown 
of the small phase error assumption. Fortunately the velocity variance is bounded 
and decreases with increasing observation time as shown next 

At low offset frequencies, the difference phase variance is given by 

(19) 

JBx 

where Bjp) are the low (high) frequency limits of integration and observation. The 
range variance is related to the difference phase variance by 
where X is the modulation wavelength. As the observation period increases, 
B\ -> 0, the range variance increases without bound. Defining Tip) = l/flip), 
the range variance may be expressed as a function of the observation time as 
ocoqt] ln(72/Ti)(A./47r)^. The velocity error is then given by 


5? ^ 2a6>;Tjln(^) 

Ti Ti 


( 20 ) 


For increasing observation time, T 2 —*■ 00 , the velocity eiTOT goes to 2 »ro. 

As a particular example, assume that an oscillator just satisfies the high fre- 
quency phase noise requirement and the phase nmsc 1//^ dependence starts at 
1 kHz. Then a < -65 dBc/Hz and (oq = 2n x I kHz. Let X = 150 mm and 
Tj = 1 ms which is appropriate for the fine ranging system at maximum satellite 
separation. Then the velocity error is given by 


or„ 




( 21 ) 


The second factor is of order unity. For example, Vln(10*°) = 4.8. Let Ti = 1 ms. 
which is a very conservative estimate of the upper frequency response fw the 
measurement Over 1 second periods, 72 = 1 s, Ae velocity error is 160/i.m/s. For 
72 = 10 s, Ae velocity error is 18/Lt/s. Over 100 s, Ae velocity enw is 2/im/s. So 
Ac l/a>^ phase noise is not likely to be a major problem. However it is also likely 
Aat Ae phase model is not very accmate over such long time periods. 


4 Conclusions 

Table 2 shows some specifications of commercially available oscillators. Hgure 2 
plots Ac phase noise spectrums. The dependence is clearly visible m Ae TRAK 
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Maker 

Frequency 

Frequency Drift 

Stability Rate 

Phase noise (dBc) at 
10 Hz 100 Hz IkHz 

TRAK 

2 GHz 

5 X 10"* 5 X 10"Vyear 

-55 


-105 

Efratom/Ball 

10 MHz 

5 X lO-'Vs 

-74 

-94 

-104 

Wenzel 

10 MHz 

1 X lO-’/day 

-89 

-109 

-119 


Table 2: Specifications of some commercially available crystal oscillators. The phase 
noise specifications are after perfect multiplication to 2 GHz. 



Figure 2: Phase noise power spectral density for three commercial osciUators: One 
from TRAK, Efratom/Ball, and Wenzel. Phase noise adjusted to 2 GHz finequency when 
required. Also shown are the fits to the lowest frequency values. 
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and Wenzel oscillators. All three oscillators easily satisfy the high frequency 
phase noise requirement Also, extending the low frequency fits to 1 kHz gives 
a < — llOdBc for all three which yields velocity errors mote than 100 times 
smaller than in the example of the previous section. 

Assuming that the frequency stability and drift rates can be extended linearly 
to different time spans, either the Efttom/Ball or the Wenzel oscillator will satisfy 
or nearly satisfy the GAMES oscillator performance requirements. However more 
information is needed on the stability to be sure. 


Derivation of stability equations 

These appendixes derives the equations used in section 2 to compute the master 
oscillator frequency and phase stability requirements. Rrst the frequency stability 
required for the velocity measurement is computed. The frequency variations 
may be either long term drift or random fluctuations. Then the phase stability 
requirement is derived. Finally, the effect of frequency variations on the absolute 
range measurement is derived. This will bound the absolute inter-satellite range 
measurement over the life of the GAMES mission. 


A Derivation of required frequency stability 

First die effect of frequency drift on the velocity measurement is determined. 
Assume that the frequency is exactly fc and the relative velocity is zero. Then 


Aor = 


Infclr 

c 


If the range changes by Ar then the new difference phase will be 


(22) 


Aa -I- Aa,. 


2nfc2[r A,] 
c 


(23) 


where A«r is the change of phase due to a range change and is given by 2nfe2Ar /c. 
If the frequency were to change instead from fe to fc-\- Afc the new difference 
phase would be 


Aa -I- Aa/ = 


2n[fc + A/,]2r 
c 


(24) 
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where Aa / is the change due to the frequency change and is given by 2n c. 
For Ar to be the resolution limit, must have Aa/ < Aa, or 


^ < — 
fc~r 


(25) 


as claimed in item 1. 


B Derivation of maximum frequency drift 


To measure the range rate to a specified accuracy requires that the error in the range 
measurement grow no more rapidly than the desired range rate. Assume that it 
is desired to measure the range rate over a time period Ar to an accuracy of Av. 
Two range measurements, each with accuracy Ar, are required. If Ar is the range 
measurement standard deviation, then the velocity enw is Av = JlArj At. If 
Ar is the peak error, then the velocity error is Av = lArjAt. Because the 
latter provides a tighter bound, it will be used. Then the range error permitted for 
each measurement will increase as Ar = ArAv/2. Assume that the oscillator 
fitequency is exactly fc- At time r, the difference phase is measured as 


Aa = 


2nfc2r 

c 


(26) 


Assume that during the time interval t to / + Ar the range increases to r + Ar. 
The new difference phase measured at time r + Ar is 


Aa + Aa,. 


2itfc2[r + AvAr] 
c 


(27) 


fw a change of Aa, = 2nfcAvAt/c and a time rate of change of the difference 
phase of 

^ = ^2Av (28) 

At c 

If the range did not change, but the oscillator frequency changes from /« to /c + A/^, 
then there will be a difference phase change to 


2rr[/, + A/J2r 

/ = 


Aa + Aa 


c 


(29) 



Therefore the time rate of change of the difference phase is given by 


Att/ _ A/- 

At cAf 


( 30 ) 


Fot the desired phase shift Aa, can be observed, it must be true that Aa//At < 
Attr/ At or 


A/c 

/cAt 


Av 


(31) 


as claimed in 2. It makes no difference if the frequency drifts gradually during 
the time interval, wonders randomly during the interval, or experiences a step 
frequency change. 


C Derivation of required phase stability 

So far, attention has focused on frequency changes. Now the attention turns to 
phase noise which has two effects. In a standard Doppler radar system, phase 
noise of the master oscillator decreases the dynamic range of the doppler system. 
When a strong reflector and a weak reflector are both in the radar field of view 
and have similar doppler shifts, the weak reflector will be hidden in the reflected 
strong signal. Because GAMES has only one target in the field of view, this is not 
a problem. 

Phase noise will also introduce a random component to the phase measurement 
Fot frequencies offset from the carrier by more than «1 kHz, the received and ref- 
erence phase are uncorrelated. The phase noise must be an insignificant component 
to the total noise budget I will specify that the phase noise power spectral density 
is about the same or smaller than the thermal noise from die PMT pre-amplifier. 
However, it should be remembered that the phase noise is prt^iortional to the signal 
power and will increase with signal power or PMT gain 

The phase noise power density (dBc/Hz) should be about the same as the 
amplifier thermal noise to signal ratio (dBsignal/Hz). Xiaoli Sun reports the signal 
power firom the PMT as 5.5 x lO"*® The noise current of a 50D resistor is given 

by i* = N FkaT fRi where NF is the noise factor, here assumed 2 times. The 
noise current is 1 .6 x 10"“ A^/Hz. Therefore the thermal noise to signal ratio is - 
65 dB/Hz which is about what the phase noise at 2 GHz must be, -65 dBc/Hz. This 
specification is easy to meet Most crystal oscillators have phase noise densities 
smaller than -80 dBc, even after fi^uency multiplication to 2 GHz. 
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Sample Rate and Number of Bits 
Required for the GLAS Waveform Digitizer 

Xiaoli Sun/JHU, December 19, 1994 

A minimum sample rate and the vertical resolution is derived for the 
waveform digitizer in the GLAS receiver. It shows that a 0.8-1 Gs/s 
6 bit A to D converter is sufficient to achieve the require ranging 
performance. 

1 . The Sample Rate of the GLAS Waveform Digitizer 

The sample rate of the GLAS waveform digitizer should be chosen 
such that the ranging error caused by the finite sample rate is 
negligible as compared with other sources of errors. The required 
ranging accuracy for GLAS is 5 cm for a roughly 30 ns wide received 
pulse (3° slope, 70m diameter laser foot print, and 5 ns wide laser 
pulses). The current design of the GLAS receiver sets the sample 
rate to 3 Gs/s, which corresponds to a range bin of 5 cm. This choice 
of the sample rate is supported by the analysis given by Gardner^ , 
which says the sample interval should be about 1 / 1 0 the rms pulse 
width. Since the laser pulse width for GLAS is 5 ns FWHM, the rms 
pulse width is 2.5ns/(2ln(2))^/2- 2.13ns. Assuming the impulse 
response of the photodetector has the same pulse width, the 
minimum rms pulse width of the signal into the waveform digitizer 
is (2.1 3^x 2.1 32)1/2=3.01 ns. Therefore, the required sample rate 
would be about 3.3 GHz according Gardner's analysis. 

On the other hand, since the input signal is band limited, the sample 
rate needs only to be twice the reciprocal of the bandwidth 
according to Nyquist Theorem. The 3dB bandwidth of the GLAS 
photodetector is about 1 50 MHz, the stop band (-30dB) frequency is 
about 400 MHz assuming a >20dB/octave roll off response. 
Therefore, the sample rate needs only to be 0.8 Gs/s. Lower sample 
rate would mean less electronics, lower power consumption, and 
possibly less signal processing time. In theory, one could completely 
reconstruct the waveform by taking the Discrete Fourier Transform 
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(DFT) of the samples and then a regular inverse Fourier transform of 
the center segment of the DFT. 

If the general pulse shape is known, even fewer samples are needed 
to recover the pulse shape. Ideally, the number of samples need only 
to be equal or greater than the number of unknown parameters of the 
received pulse shape. For example, three samples would be necessary 
if the pulse shape were Gaussian with unknown time delay, width, 
and amplitude and no noise. More samples are required in practice in 
order to average out the noises. The parameter values of the 
received pulse can be obtained through a standard data fitting 
technique, such as maximum likelihood or minimum mean square 
error data fitting. In general, the timing accuracy can always be 
much smaller than the sampling interval unless the pulses are 
perfect rectangle with zero rise and fall times. As an example, the 
LeCroy LW41 0/420 arbitrary waveform generator can form a time 
shift in Its output waveform by as little as 1 00 ps, a small fraction 
of the 2.5 ns sampling interval, with 8 bit vertical resolution^. 

For GLAS, no assumption can be made about the exact pulse shape 
since they represent the features of the ground target to be studied. 
The sample rate must be determined based on the transmitted laser 
pulse shape and the receiver impulse response. Intuitively, the 
sample rate has to cover the rate of change of the received signal. 
Based on this principle, the sampling interval has to be a fraction of 
the pulse rise and fall times. Since the GLAS detector has a 3 dB 
bandwidth of 1 50 MHz, the rise time of the output signal Is at most 
0.35 times the reciprocal of the bandwidth, or 2.3ns. We can safely 
set the sample rate to 850 Ms/s which gives two samples per rise 
time. This sample rate is also close to that given by the Nyquist 
Theorem. One could reconstruct the received pulse to a great 
accuracy by connecting the samples using standard linear or cubic 
spline interpolation. The centroid of the received pulse which gives 
the range of the ground target may be calculated approximately by 
summing all the samples weighted by their arrival times. One may 
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further improve the ranging accuracy by using the Simpson's rule in 
the summation, as in the theory of numerical integration. 

In conclusion, the 3Gs/s sample rate for the GLAS waveform 
digitizer is shown to be excessive and the suggested is 0.8-1 Gs/s. 
The analysis given by Gardner about the quantization error appears 
to be inaccurate. An exact analysis of the ranging error due to the 
sampling rate seems difficult. A numerical solution may be obtained 
through a computer simulation using the GLAS simulator. The 
ultimate solution has to come from the test results of the GLAS 
breadboard receiver. 


2. The Number of Bits of the GLAS Waveform Digitizer 

The finite step size of the A/D converter causes a quantization error 
In each sample of the output. The quantization errors may be 
modeled as identically and independently distributed random 
variable and they all follow a uniform distribution of the digitizer 
step slze3. 

The output of the waveform digitizer may be written as 
y(n) = x(nAi) + = x(n) + e- 

with x(t) the input analog signal. At the sampling Interval, and ej the 
quantization error. The measured target range is proportional to 
centroid of the received pulse, which can be written as 

^ jA/y(/) ^ i x{i) + ^ i e, 

i i I 
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Assuming the quantization errors are small compared to the average 
pulse amplitude, i.e. , e, x(/) « i , the pulse centroid may be 

t / i 

approximate as 
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The quantization error normalized with respect to the rms pulse 
width, Op, is given by 


Vvar{r,} ^ 1 
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where Ax=Xmax/2N is the step size, Xmax is the full scale voltage, and 
N is the number of bits. 


We further assume that the received pulse shape is Gaussian, i.e., 
x(n)=xoexp{-(nAt)2/2op2} with xq the peak pulse amplitude. The 
summation In the centroid calculation is assumed to extends to ±2op. 
The scaling factor of the waveform digitizer is chosen such that 
Xmax=axo with a >1 the dynamic range. Under these conditions, the 
quantization error becomes 
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To calculate the required number of bits for GLAS, we assume 
Op=l 3ns (30ns FWHM pulse width) and At=1 ns (IGs/s). The 

quantization error for a 8 bit A/D converter with 3dB dynamic range 
(a=2) is e»5.76xl 0“'^. 
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The relative ranging error required by GLAS is 


2 ^ 


The quantization error has to be much smaller than emax- If we 
choose e=emax/l 0, the number of bits required is 6 bits. An 8 bit A/D 
converter will give a dynamic range of a =8 (1 8dB). 

The desired dynamic range for GLAS receiver is 26dB, which 
accounts for a factor of 1 0 variation in the ground reflectivity and a 
factor of two APD gain fluctuation. Therefore, an automatic gain 
control (AGC) amplifier is still required in front of the waveform 
digitizer. 
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